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ABSTRACT 

Doubly-degenerate binary systems consisting of two white dwarfs both composed of 
carbon and oxygen and close enough that mass is transferred from the less massive to 
the more massive are possible progenitors of type la supernovae. If the mass transfer 
rate is slow enough that the accreting white dwarf can reach a mass of 1.38 M0 then it 
can ignite carbon degenerately at its centre. This can lead to a thermonuclear runaway 
and thence a supernova explosion. However if the accretion rate is too high the outer 
layers of the white dwarf heat up too much and carbon ignites there non-degenerately. 
A series of mild carbon flashes can then propagate inwards and convert the carbon 
to neon relatively gently. There is no thermonuclear runaway and no supernova. We 
examine the critical accretion rate at which ignition switches from the centre to the 
surface for a variety of white dwarfs and find it to be about two fifths of the Eddington 
rate. In a real binary star the mass transfer rate falls off as mass transfer proceeds and 
the system widens. Even if the initial transfer rate is high enough for carbon to ignite at 
the outer edge, if such rapid accretion were to persist, we find that it can extinguish if 
the rate drops sufficiently quickly. The interior of the white dwarf remains carbon rich 
and, if sufficient mass can still be transferred from the companion, it can eventually 
ignite degenerately at the centre. The primary white dwarf must be about 1.1 Mq or 
above and the companion about 0.3 M©. Though white dwarfs of such low mass are 
expected to be pure helium we note that a star of initial mass 2.5 M© has a CO core 
of about 0.3 Mq when it begins to ascend the asymptotic giant branch. Alternatively 
if the accretion rate can be limited to a maximum of 0.46 of the Eddington rate then 
a 1.1 Mq white dwarf accretes sufficiently slowly to explode from a companion white 
dwarf of any large enough mass. 

Key words: white dwarfs, supernovae: general, stars: evolution, binaries: close. 



1 INTRODUCTION 

Type la supernovae (SNe la) are the brightest objects in nor- 
mal galaxies. They appear to be standard candles an d hence 
are u seful cosmological measures of the Universe jBaadd 
[W3^. They are also a major source of iron iTout et"aLl20Ml) 
which tends to be trapped in the neutron star remnants 
of other supernovae. Ho wever^ their progenitors remain un- 
certain (see for example iToutI l|200^) and this leaves their 
standard nature questionable particularly in the light of the 
chemical evolution of the Universe. 

Type I supernovae have no hydrogen lines and type la 
are further distinguished by their prominent silicon lines. 
They are almost certainly exploding carbon and oxygen 
(CO) white dwarfs. Their available nuclear energy exceeds 
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their binding energy so that the whole star can be de- 
stroyed in a thermonuclear runaway. Most of the material 
approaches nuclea r statistical equili brium and about 0.6 M© 
of '^'^Ni is expelled iRopke fc Hillebr andt 20q1). The radioac- 
tive decay of ^^Ni to ^^Fe via ^^Co powers the su pernova 
and the decay signature has been clearly observed iBranchl 
1998). It is fairly certain that SNe la are detonated by mass 
accretion on to the CO white dwarf. As the mass of a cold 
white dwarf increases towards the Chandrasekhar mass, the 
gravitational collapse heats the material and fusion ignites 
in the degenerate core. Typically CO white dwarfs explode 
at about 1.38 Mq. 

Accreting white dwarfs have been known for some time 
as the engines of cat aclysmic varia bles, the source of no- 
vae and dwarf novae llWarneilll995l) . and so were the first 
candidates to be considered. However if the accreting ma- 
terial is hydrogen-rich, accumulation of a layer of only 
10~^ — 10"'^ Mq of cold material leads to degenerate igni- 
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tion of hydrogen burning sufficiently violent to eject most, 
if not all of or more than, the accreted layer in the well 
known nova outbursts of cataclysmic variables. The white 
dwarf mass does not significantly increase and ignition of 
its interior is usually avoided. However if the accretion rate 
is high, M > lO"'^ Mq yr""*^, compressional heating of the 
surface layers raises the degeneracy and hydrogen can burn 
relati vely gently as it is accreted, bypassing novae explo- 
sions l|Paczvhski fc Zvtkow , 19781. allowing the white dwarf 
mass to grow. Though, if it is not much larger than this, 
M > 3 X 10"'^ Mq yr~^, hydrogen cannot burn fast enough 
so that accreted material builds up a giant-like envelope 
around the core and burning shell which eventually leads 
to more drastic interaction with the companion and prob- 
ably the end of the mass transfer episode. Rates in the 
narrow range for steady burning are found only when the 
companion is in the short-lived phase of thermal-timescale 
expansion as it evolves from the end of the main sequence 
to the base of the giant branch. Super-soft X-ray sources 
iKahabka fc van den Heuvellll997l) are probably in such a 
state but cannot be expected to remain in it for very long 
(a SNe la rate of only about one millionth of the observed 
rate form by this channel in the population synthesis of 
iHurlev. Tout fc Polj (|20o3)) and white dwarf masses almost 
never increase sufficiently to explode as SNe la. If He-rich 
material is accreted instead about 0.1 — 0.15 Mp, of degener - 
ate material can accumulate before ignition llNomotdll982ll 
and an accretio n rate above 3 x Mp^yT~ ^ can lead to 
steady burning llKawai. Saio fc Nomotdligsil . 

The currently popular model overcomes these problems 
by postulating that, when the mass-transfer rate exceeds 
that allowed for steady burning, only just the right fraction 
of the mass transferred is actually accreted by the white 
dwarf. A viable mechanism for this is a strong wind from the 
accretion disc that expels material from the system befor e 
it reaches the white dwarf iHachisu. Kato fc Nomotolll99d) . 
Alternatively the white dwarf might indeed swell up to gi- 
ant dimensions but the resulting common-envelope evolution 
could be sufficiently efficient that the small amount of excess 
material can be ejected without the cores spiralling in. This 
is qui te c onsistent with the find ings of iNelemans fc Tou3 
JioO^ and lNelemans et al.l ll2000l) that such efficiency neces- 
sary for at least one phase of common-envelope evolution in 
the formation of close double white dwarf systems. Helium- 
accreting CO white dwarfs were once popular when it was 
thought that a thermonuclear runaway in the accreted he- 
lium layer could set off the CO core in an edge-lit detona- 
tion (IWooslev fc Weave jlT993) . However 2-D models indi- 
cate that central ignition is unlikely and light curves and 
spectra do not fi t well so they are n o longer considered vi- 
able progenitors iBranch et alJ[T995^ . 

In this work we consider the alternative that many 
complications can be avoided altogether if the white dwarf 
can accrete material of similar composition to itself, car- 
bon and oxygen. This can be achieved in a doubl e degener- 
ate m odel in which two CO white dwarfs merge JWebbinkI 
Il984} ). Both form as the cores of asymptotic giant branch 
stars and, after one or two phases of mass transfer, one or 
both of which involves substantial orbital shrinkage through 
common-envelope evolution, they are brought close enough 
together that angular momentum losses owing to gravita- 
tional radiation lead to Roche lobe overflow. A major diffi- 



culty with this model lies in the fact that high accretion rates 
also heat the white dwarf. This heating is due mainly to the 
gravitational compression of the white-dwarf material as the 
mass increases. The actual accretion luminosity liberated by 
the material falling down the potential well on to the white- 
dwarf surface has a negligible heating effect in comparison 
and is in any case mostly liberated in an accretion disc or 
boundary layer and radiated away. A white dwarf can lose 
heat by conduction to the surface and radiation or by in- 
terior neutrino loss processes such as electron-positron pair 
production and annihilation. The material at the surface un- 
dergoes the most rapid compression as it becomes degener- 
ate. At slow rates energy released in these surface layers can 
escape fast enough but at high rates the surface layers reach 
temperatures required for fusion to begin. True ignition oc- 
curs when energy production by carbon fusion exceeds that 
lost in neutrinos, a situation that is rapidly followed by the 
onset of convection but at velocities that are not sufficient 
to carry aw ay igniting packets bef ore thermonuclear run- 
away ensues. iNomoto fc IbenI lligSST) calculated that carbon 
ignites near the surface at constant accretion rates in ex- 
cess of about one fifth of the Eddington rate (see section [IJ. 
Ignition at the surface raises the degeneracy so that burn- 
ing of carbon to neon and sodium proceeds relatively gently. 
Surface burning triggers the gentle ignition of a deeper shell 
and carbon is thus burnt successively throughout the white 
dwarf. Once the central carbon has been exhausted the core 
can't ignite until hot enough for neon burning by photo- 
disintegration or oxygen fusion but by this point electron 
capture by magnesium has taken the collapse beyond the 
point at which the available nuclear energy can explode the 
white dwarf. It simply collapses quietly to a neutron star 
releasing energy in neutrinos. Thus for a SN la a thermonu- 
clear runaway must begin in the core before carbon ignites 
sufficiently at the surface. 

Because white dwarfs expand as they lose mass, unless 
the Roche lobe of the mass donating companion expands 
even faster the very process of mass transfer causes the mass- 
losing (lower mass) white dwarf to overfill its Roche lobe yet 
more and mass transfer accelerates to dynamical rates. Only 
when the mass ratio q = Mioscr/A/accrotor < 0.628 can this 
positive feedback be avoided. It is not known what the final 
outcome of dynamical mass transfer is but it is likely that 
much of the loser is lost to the interstellar medium and that 
the remainder is accreted at very high rates. This already 
excludes the major ity (about four-fifths in th e population 
synthesis models of lHurlev. Tout fc Pols! l|20o3)) of mergmg 
CO white dwarf pairs which tend to have component masses 
between 0.5 and I.IM©. In the remaining cases, where dy- 
namical mass transfer is avoide d, the initial rate is still in 
excess of iNomoto fc Ibeiil lll985ll 's limit and so merging CO 
white dwarfs are not currently favoured SNe la progenitors. 

Various me chanisms might force merging white dwarfs 
to spin rapidly. iPiersanti et alJ ll2003allbl) modelled the ef- 
fects of rapid spin on carbon ignition in CO white dwarfs. 
They found that compressional heating was not much af- 
fected but that the rate of diffusion of heat to the centre was 
slowed. This generally reduces the maximum rate of accre- 
tion that allows central carbon ignition. On the other hand 
spinning near breakup prevents accretion but also distorts 
the white dwarf so that it spins down by its own gravita- 
tional wave emission, while accreting just enough mass to 
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keep it spinning close to break up. This leads to a natu- 
ral limit to the accretion rate of about 4 x lO~^M0yr~^, 
low enough to allow central carbon ignition. In this work we 
only consider slowly rotating white dwarfs because we are 
interested in dynamically stable mass transfer. For instance, 
when a 0.3 Mq white dwarf begins mass transfer to a 1.1 Mq 
white dwarf (as in section|SJ the period of the binary is 146 s, 
well above the break-up spin of the 1.1 Mq white dwarf of 
5.6 s. It might however be argued that accretion from the in- 
ner edge of a Keplerian accretion disc would still spin up the 
massive white dwarf before it reached its ignition mass. How- 
ever this is avoided if the white dwarf has a magnetic field 
strong enough to disrupt the accretion disc as in the pola r 
and intermediate-polar cataclysmic variables JWarneill995l) . 
Indeed more massive wh ite dwarfs are thought to form wit h 
stronger magnetic fields iWickramasinghe fc Ferrariol20o3) . 
We thus expect the accreting white dwarfs which interest us 
here to be relatively slowly rotating and so can be modelled 
without rotation. 

Here we reconsider the limit on the accretion rate. After 
describing the details of our evolution code and how wc iden- 
tify carbon ignition (section I^J we first wo rk wi th cons tant 
rates (sectional for direct comparison with lNomoto fc IbenI 
il985() . With modern computer power we can calculate on 
a much finer grid of initial masses and accretion rates. Sec- 
ondly, we consider accretion at a constant fraction of the 
Eddi ngton rate (sectional and find somewhat higher limits 
than lNomoto fc IbenI ||l983) . Finally we consider true binary 
evolution in which the rate is initially high but falls off as the 
donor white dwarf loses mass and the system widens (sec- 
tion 01 • In section |S| we explore the range of initial binary 
star parameters for which we find central carbon ignition 
and thence a supernovae and in section Q we present our 
conclusions and consider whether realistic progenitor sys- 
tems might exist. 



2 STELLAR EVOLUTION MODELS 

We use the Cambridge STARS code to make, evolve and 
accrete on to white dwarfs. The STARS code is th e most re- 
cent version of th e Eggleton evolution program jEggletonl 
Il97ll. Il97l I1973D . The equation of state, which includes 
molecular hydrogen, press ure ionizat i on an d Coulomb in- 
teractions is discussed bv iPols et al.l 1^9^. We take the 
initial composition of the star to be uniform with a hy- 
drogen abundance X = 0.7, helium Y = 0.28 and metals 
Z = 0.02. The metal mixt ure is according to the met eoritic 
abundances determined bv lAnders fc Grevessd ^1989^ . Only 
nuclear burning that affects the structure is included, hydro- 
gen burning by the p-p chain and the CNO cycles, helium 
burning by the triple-a reaction and reactions with ^'^C, ^^N 
and ^''O and carbon burning via 

12^ + only. Other ISO- 
topes and reactions a re not explicitly modelled. T he reaction 
rates are taken from ICaughlatL-fcjQwleJ <ll98j) and opac- 
ity tables are from llglesias. Rogers fc WilsonI il992l) and 
[Alex ander fc Ferguson' (1994). An Eddington approxima- 
tion ( W ooUey & Stibbs 1953) is used for the surface bound- 
ary conditions at an optical depth of r = 2/3. 

We initially construct white dwarf models from those of 
asymptotic giant branch stars with appropriate core masses 
by stripping them of their envelopes at very high mass-loss 
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Figure 1. The internal temperature density profile of a white 
dwarf of initial mass 1.2 Mq with a central temperature of lO'^'^ K 
made from the core of a 5Mq star. The solid line is this initial 
state before accretion. Nearly all the mass lies in the isothermal 
core. Most of the variation in temperature occurs through the 
very thin atmosphere. The dotted line has accreted to 1.38 Mq 
at 2 X 10~^ Mq yr~i and is about to ignite carbon at its centre. 
The dashed line has accreted to 1.35 Mq at 3 x 10~^MQyr~l 
and is about to ignite carbon just below its surface. The ignition 
of each of these is shown in more detail in figures 1^ and O 

rates to mimic the binary formation process. This leaves CO 
cores which can cool, without further mass loss, to become 
white dwarfs. We obtain white dwarfs of different temper- 
atures by varying the cooling time and we can form white 
dwarfs with different compositions by evolving to the AGB 
from main-sequence stars of different initial mass (see tables 
□ andHll. 

Once we have a white dwarf we can accrete mass back on 
to it. However the mass we accrete now is of the same com- 
position as the surface of the white dwarf. There are then no 
problems with numerical diffusive mixing. Because our evo- 
lution code is hydrostatic, to ensure numerical convergence, 
we must prime our white dwarfs for rapid accretion by in- 
creasing the rate from zero to the desired amount in several 
short steps. These steps are short enough that this is effec- 
tively instantaneously compared with the thermal timescale 
in the regions of interest. 

When carbon ignites in a thermonuclear runaway it 
manifests itself by a sudden rise in carbon burning luminos- 
ity followed by breakdown of the evolution code. Figure 
shows how the internal run of temperature and density from 
the cool surface to hot relatively isothermal interior changes 
during accretion at rates just above and just below the crit- 
ical rate at which ignition moves from the centre to the 
surface. To be sure of what has happened, we determine a 
carbon ignition curve as a function of temperature and den- 
sity by equating the carbon burning energy production with 
that lost in neutrinos - see Figures |2| and |2l Above and to 
the right of this line we have the right conditions for carbon 
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Figure 2. The line crossing both axes in the carbon ignition curve 
on which energy generation by carbon burning equals the energy 
loss in neutrinos. The other lines are the run of temperature and 
density through the white dwarf, at late times, increasingly later 
from bottom to top, beginning with the dotted line of figure^ In 
this case these cross the ignition curve at the centre and a ther- 
monuclear runaway can begin there under degenerate conditions. 
The accretion rate is 2 X 10"'' MQyr"'^. 
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Figure 3. As figure |21but be ginning with the dashed line of 
figure Q Ignition occurs close to the surface and carbon ignites 
under much less degenerate conditions. The accretion rate is 3 X 
lO~^M0yr~^. Grey filled circles indicate the path of the mass 
shell where carbon first ignites at mass 1.354 just 0.008 Mq 
below the surface. 



ignition to run away. Below it, any carbon burning ceases 
if the accretion is turned off. On the same axes in figures 2 
and 3 we plot the internal temperature against density for 
the models leading up to ignition. If a line from our inter- 
nal model crosses the ignition curve then we know that the 
carbon in the star has ignited at that particular point in the 
model. There is then a clear distinction between those mod- 
els that ignite at the centre (figure|5J and those which ignite 
in the outer layers just below the surface (figureEJ. There is 
not a smooth transition of the ignition point from the centre 
to the surface because, in both cases, the intervening parts 
of the white dwarf lie well away from the ignition curve. 
The transition is sudden and so the critical rate at which it 
occurs is well determined for each of the cases we consider. 
We note further that this ignition in the outer layers differs 
from the concept of an off-centre but still relatively central 
ignition in the degenerate interior which could still lead to 
thermonuclear runaway and supernovae. 



3 CONSTANT ACCRETION RATES 

iNomoto fc IbenI il985l) considered accretion at constant 
rates. For a direct comparison we first examine a range of 
initial conditions and constant accretion rates (as in figures|5| 
and |2|to estimate the critical rate at which ignition switches 
from centre to surface and how it depends on initial temper- 
ature and central carbon to oxygen ratio. Table^shows the 
critical accretion rates at which the ignition switches from 
the centre to the outside of the star for different white dwarf 
models. We see that, in general, the hotter the initial central 
temperature, the higher the critical accretion rate. The hot- 
ter central temperature makes it easier for carbon to ignite 
at the centre before the outside. For white dwarfs of similar 
temperature, the higher the initial mass, the lower the crit- 
ical accretion rate. Also, the higher the carbon to oxygen 
ratio, the lower the critical accretion rate. 



4 EDDINGTON ACCRETION RATES 

Because it has some physical meaning and could place a real 
limit on the accretion rate, it is also interesting to investigate 
the response of a white dwarf to accretion that varies at a 
fixed fraction of the Eddington rate MRnn. lNomoto fc IbenI 
found the critical rate to be about one fifth of Medd. 
In the case of one dimensional spherical accretion the kinetic 
energy of the accreting material is liberated in a shock at 
the surface of the white dwarf and, if radiative transfer is 
the only mechanism for the energy to escape, the consequent 
radiation pressure slows the infalling material. A maximum 
accretion rate, Medd, is reached when the radiation force on 
particles balances gravity at the surface of the white dwarf 
so 

Medd = Ri, (1) 

where Ri is the radius of the accreting white dwarf, c is the 
speed of light and k is the opacity of the accreting material 
which is dominated by electron scattering. This is found by 
equating the Eddington luminosity and the accretion lumi- 
nosity, that released by material falling from infinity to the 
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Table 1. Critical constant accretion rates. Column 1 is the original zero-age main-sequence mass of the star from which the white dwarf 
was made. Column 2 is the white dwarf mass. It's central carbon and oxygen abundances are given in columns 3 and 4. The initial 
central temperature Tc is in column 5. Column 6 is Merit the critical accretion rate at which the ignition switches from the centre to the 
outside of the white dwarf. 



surface of the white dwarf. The Eddington limit decreases 
as Ml increases because the radius R\ falls as Mi rises. 

In all cases accretion at Medd is sufficiently fast to 
cause ignition in the outer layers and so we also consider 
accretion at a fraction of A^edd. In spherically symmetric 
accretion Medd represents t he absolute limit to the rate at 
which mass can be accreted, tean fc Webbinkll999lh make a 
more careful analysis, taking into account the potential dif- 
ference between the inner Lagrangian point and the white 
dwarf surface and so permit higher accretion rates. If ac- 
cretion is not spherically symmetric so that material can 
accrete along some directions and radiation escape in others 
even higher accretion rates could be possible. On the other 
hand, if the opacity k increases rapidly as material is driven 
off, accretion may be limited to a lower rate in a similar 
manner to the limit on hydrogen accretion in the model of 
[kachisu. Kato & Nomoto (1996). 

The radius Ri_ depends on the mass Mi of the ac c reting 
white dwarf according to the formula of iNauenberd (Il972l) 
as 

where M^h = 1.44 M© is the Chandrasekhar mass. For a 
selection of white dwarfs table |5| lists the critical fraction 
of Medd at which we move from central to outside ignition 



and Figure 3] shows how this critical rate varies with time 
for each case. We see that we can accrete up to about two 
fifths of Medd and stil l find central ign i tion. T his is double 
the rate estimated bv iNomoto fc IbenI il985l) . The higher 
the initial mass of the white dwarf the higher the critical 
fraction of the Eddington rate. We do not profess to have 
identified an actual physical mechanism that would limit the 
accretion in this way but rather wish to indicate what would 
be necessary. 



5 ROCHE LOBE OVERFLOW IN WHITE 
DWARF SYSTEMS 

The actual rate of mass transfer between double white 
dwarfs often exceeds A/edd but it is not constant. Rather 
it falls off as mass is transferred and the system widens in 
response to the expansion of the radius of the mass-losing 
(donor) white dwarf. The separation is controlled by two 
competing processes. Gravitational radiation acts to shrink 
the orbit while mass transfer from the less to more massive 
star acts to widen it. For stable mass transfer the Roche lobe 
must expand faster than the white dwarf, which grows in re- 
sponse to mass loss. The Roche lobe itself grows because the 
system expands as mass in transferred to the more massive 
companion. The combination of these changes determines 
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Table 2. Headings as for table Q except the last column whicli is tlic fraction of the Eddington accretion rate for which the carbon first 
ignites in the outside of the white dwarf. 
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Figure 4. The variation with time of accretion rate at the critical 
fractions of the Eddington rate for the four cases listed in tablel^ 
The rate is proportional to the radius of the white dwarf and so 
falls off as mass increases. The top line is for accretion on to the 
white dwarf of initial mass 1 Mq and the bottom of 1.3 Mq 



the mass-transfer rate necessary to maintain R2 = Rl2 and 
the so the system widens. As it widens the system experi- 
ences weaker gravitational-radiation braking and the mass- 
transfer rate falls. 

We set up a simple model that allows us to determine 
this mass-transfer rate and limit it to a maximum fraction 
/ of Medd. In this way we can find systems in which the 
primary white dwarf ignites at the centre even though the 
initial mass transfer rate might be considerably higher than 
the critical rates calculated earlier. 

In a binary system the Roche lobe radius i?Li is the ra- 
dius of a sphere which encloses the same volume as that en- 
closed by the last sta ble equipotentia l surface around star i 
(i = 1, 2) as given by jEggletonlllOSll . 

^ 0.49gf/^ 

where qi is the mass ratio Mi/M^-i and a is the separation 
of the two stars. This formula is accurate to 1 percent for 
all q. 

In a binary system of two white dwarfs, when the less 
massive, of mass M2 and radius R2, overfills its Roche lobe 



it loses mass, some of which is accreted by its companion, 
of mass Ml and radius R\. For stable mass transfer the 
rate (—M2) is determined by the angular momentum loss 
from the orbit because M2 adjusts to keep R2 ~ Rh2 and 
R2 ~ Rh2- We model the mass-loss rate from the donor 
white dwarf by 

M2 = -Moe^, (4) 

where AT? is the amount by which the white dwarf overfills 
its Roche lobe, 

AR = R2- Rh2 (5) 

and H is the p ressure scale height (|RrtteJ ll98Sl : 
iMartin fc Toull2005l) . This formula takes account of the fact 
that stars have a thin atmosphere above their photosphere 
and Roche lobe overfiow begins while the photosphere it- 
self is still below the inner Lagrangian point. It also allows 
us to make a smooth transition to stable mass transfer as 
the white dwarfs are driven together. Because it is a steep 
function of AR we have R2 ~ Rh2 once mass transfer has 
begun. The constant Mo is the mass-transfer rate when the 
donor star fills its Roche lobe exactly so that AT? = 0. We 
use Mo = 3 X 10~^. This is approximately the critical mass 
transfer rate where ignition moves from the centre to the 
outside of the WD. If this rate is exceeded then we slightly 
over estimate the mass transfer rate because the stars are 
slightly closer than they would be if AT? = and |Jgr| (see 
equation IHl below) is slightly larger. If M < Mo we slightly 
underestimate the rate. Figure |^ below shows that this rate 
is typical of the systems in which we are interested. 

We take TT to be a constant fraction of one thousandth 
of the radius of the losing white dwarf. Such a small fraction 
is typical of white dwarfs so that T?2 is indeed very close to 
T?L2. Any variation in H is unimportant because we are only 
interested here in the equilibrium mass-transfer rate. 

We limit the mass accretion rate to /Medd, where < 
/ < CO is a constant. A factor / 2> 1 allows all the mass 
transferred to be accreted (ATi = —M2) while / = 1 applies 
the Eddington limit. If the mass-transfer rate is less than 
fMEUu all of the mass lost by the donor is accreted on to 
the companion. In general the accretion rate on to star 1 is 
given by 

Ml = min(/MEDD, -M2). (6) 

When the limit is active there is mass loss from the system 
because star 2 is losing mass more quickly than star 1 is 
gaining it and M = Mi -I- M2 ^ 0. This mass is blown away 
from the binary in a wind which carries angular momentum 
with it. We assume it has the specific angular momentum of 
the accreting white dwarf whence it is blown off and so the 
rate of loss of angular momentum is 
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Figure 5. Accretion rates on to a 1.2 Mq white dwarf. The solid 
line is from a white dwarf companion of mass 0.3 and the 
dotted line of mass 0.9 Mq. Both are limited by Medd a^nd this 
limit, which depends only on the mass of the gainer, is effective 
over the initial straighter section. Both tracks end when the accre- 
tor would have reached 1.38Mq. Note that no account of carbon 
burning has yet been taken. 




Figure 6. Outcome for a 1.1 Mq white dwarf accreting from a 
companion white dwarf of initial mass M2 with an upper limit to 
the accretion rate of /A/edd. In the single shaded area the white 
dwarf ignites at the centre. The blank area it accretes too fast and 
ignite at the outside. In the cross-hatched area the initial total 
mass M = Mi + M2 ^ 1.38 Mq so that central ignition cannot 
be reached. 



Jwind = Af (^a) n, (7) 

where il is the orbital angular velocity. For these close sys- 
tems loss of angular momentum owing to gravitational radi- 
ation is the mechanism driving the evolution. For two point 
masses in a circular orb it the rate of change of an gular mo- 
mentum jgr is given bv lLandau fc Lifshitd il95l|) 

jg, _ 32G^ MiM2(Mi + M2) 

We assume that there are no other angular momentum losses 
from the system so that the total angular momentum loss is 

We can solve the three differential equations Q, Q 
and with a Runge Kutta method and hence find the 
accretion rate. Figure JSJ shows an example of accretion 
rates on to the same mass of accreting white dwarf from 
two different mass companions. The rates both begin at the 
Eddington rate and fall off in such a way that the accretion 
rate from a lower mass companion is always less than or 
equal to that from a higher and the time spent at rates 
above the critical rates we reported in sections |3 and 2] is 
significantly reduced for the lower mass companion. They 
both switch to accreting all of the mass that is lost from the 
companion once this rate is less than the Eddington rate. 

We note here that, had we not used an Eddington lim- 
ited accretion rate, the binary of 1.2 Mq and 0.9 Mq white 
dwarfs would have begun dynamically unstable mass trans- 
fer. The limit on accretion rate has a stabilizing effect on 
the mass transfer. 



6 CARBON IGNITION IN BINARY STAR 
MODELS 

We now consider detailed evolution of the accretion at vari- 
able rates according to the binary model. We begin with a 
separation such that the lower-mass white dwarf is about 
to fill its Roche lobe. The initial white dwarfs are those de- 
scribed in section Eland used in section [l] In order for accre- 
tion to switch on smoothly by equation ||1J we must prime 
our white dwarf model by rapidly building up the accretion 
rate over several short time steps to about 10~^MQyr~^ 
as described in section |5| This means our surface layers are 
slightly hotter than they should be at the onset of accretion 
but the effect is soon swamped and could in any case only 
reduce the likelihood of central ignition. 

Figures |S| and [7| illustrate these, our most realistic mod- 
els, for initial mass accretors of 1.1 and 1.2 Mq of the com- 
positions given in table |21 If the total mass of the system is 
less than 1.38 Mq central ignition can't occur so the cross- 
hatched region is always ruled out. To the right of this the 
models in the single shaded regions show central carbon ig- 
nition while those in the unshaded regions ignite at the sur- 
face. For a white dwarf of mass 1.1 Mq we find central igni- 
tion for any mass companion with an accretion rate 46 per 
cent of the Eddington limit. Similarly for a 1.2 Mq accreting 
white dwarf we can succeed with any mass companion with 
accretion limited to 40 per cent of the Eddington rate. The 
curves rise steeply as / approaches 1 and the mass trans- 
fer can be fully conservative for lower-mass companions and 
still lead to central ignition. 

Figure |H| shows the highest mass of the donating white 
dwarf for which we find central ignition for a given mass of 
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Figure 7. Labels as for figure 151 for a 1.2 Mq accreting white 
dwarf. 




Figure 8. The maximum mass of donor white dwarf M2 for which 
we find central ignition of a white dwarf of initial mass Mi with 
fully conservative mass transfer / 3> 1- Circles are the masses 
actually tested. Systems in the single shaded region can reach 
central ignition. Those in the cross-hatched region do not have 
enough mass. Note that the result is very similar if / = 1. 



the accreting white dwarf with fully conservative evolution 
(/ ^ 1). Above this mass the mass-transfer rate is too high 
and we find carbon ignition at the outside and consequently 
no supernova. 



7 CONCLUSIONS 

Using detailed models of accreting CO white dwarfs we 
have critically investigated what rates allow central degener- 
ate carbon ignition as opposed to ignition near the surface. 
To compare with earlier work we first considered accretion 
at constant rates and then at fractions of the Eddington- 
limited ra te. Ou r resul ts are qualitatively similar to those of 
^omoto fc IbenI 1^8^ but we find central ignition at rates 
up to twice what they did or two-fifths of A/edd- This can be 
attributed to their relatively course grid and updated stellar 
physics. We have also considered real double white dwarf bi- 
naries with mass transfer by Roche lobe overflow. Variation 
in the initial temperature and composition of the accreting 
white dwarf have only a small effect. We note that with fully 
conservative mass transfer a 1.1 Mq CO white dwarf with a 
companion CO white dwarf of 0.3 A/q can reach the condi- 
tions required for a SN la even though the accretion rate ini- 
tially rises to a maximum of almost 1.7 x 10~^ MQyr~^ and 
carbon burning in the non-degenerate layers reaches a lumi- 
nosity of nearly 21 Lq. The mass transfer rate and hence the 
luminosity fall off sufficiently quickly with the mass-transfer 
rate that the burning is extinguished before it runs away 
and burns the entire white dwarf. 

Though standard single star evolution with reason- 
able mass loss only produces CO white dwarfs from about 
0.55 Mq to 1.2 Mq we note that a star of initial mass 2.5 A4q 
has a CO core of 0.28 Mq at the start of the asymptotic giant 
branch, immediately after convective core helium burning. 
We also recall that one or two phases of common envelope 
evolution are required to bring the two white dwarfs close 
enough for gravitational radiation to operate and that it 
would be easy to strip such a star's hydrogen rich enve- 
lope to leave a naked helium star of about 0.55 Mq. Such 
a star is very small and would not undergo further interac- 
tion until gravitational radiation sets in but it is also very 
luminous (about 200 Lq) and maintains this luminosity for 
a long time (about 5 x 10* yr owing to helium burning. This 
could well drive off the helium rich material, in a similar 
way to Wolf-Rayet stars, and expose the required l ow-mass 
CO core. lHamann. K oestcrk e fc Wessolowskil (Il995l) suggest 
a mass-loss rate from naked helium stars of as much as 
^Q-ii.95(^/^^-)i.5 ^jjjgjj jg considerably more than would 

be needed here. 

Undoubtedly the parameter space leading to these 
conditions is smal l . In the population synthesis of 
iHurlev. Tout fc Pols! ll2002ll only one in flve hundred of the 
merging CO white dwarfs would actually make it to cen- 
tral carbon ignition under these conditions and thence a 
supernova rate of only about a thousandth of that observed. 
However much of the physics, including that of common en- 
velope evolution and mass loss from naked helium stars is 
highly uncertain, as are the initial mass functions, mass ra- 
tio distributions and separation distributions of binary stars 
so any population synthesis model is still very uncertain. We 
note also that interactions between stars in dense stellar en- 
vironments can significantly increase the number of merging 
white dwarfs. Jhara & Hurley ( 200i) find a factor of fifteen 
increase in open clusters. 

On the other hand if the mass accretion rate can be 
limited to 0.46 of M^dd then a 1.1 Mq white dwarf accreting 
from a companion of any mass above 0.28 Mq can make 
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it to a supernova. This would in deed be the case if a disc 
wind s imilar to that postulated bv lHachisu. Kato fc Nomotol 
l)l99(tl operated in the double-degenerate case in a similar 
way to the single-degenerate model. 

We conclude that double CO white dwarfs remain as vi- 
able progenitors of SNe la as any others currently proposed 
but are unlikely to be responsible for the majority of SNe la. 
In particular, if the accretion luminosity were somehow lim- 
ited to less than 46 per cent of the Eddington rate in the 
early stages, then stable mass transfer could account for an 
important fraction of observed Sne la. 
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